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Abstract 

The Bg differential production cross section is measured as functions of the trans- 
verse momentum and rapidity in pp collisions at = 7 TeV, using the Bg — > ]/tp (p 
decay, and compared with predictions based on perturbative QCD calculations at 
next-to-leading order. The data sample, collected by the CMS experiment at the LHC, 
corresponds to an integrated luminosity of 40 pb _1 . The Bg is reconstructed from the 
decays ]/tp — > and (p — > K + K~. The integrated Bg cross section times Bg — > ]/ip(p 

branching fraction in the range 8 < pj < 50 GeV/ c and |y B | < 2.4 is measured to be 
6.9 ± 0.6 ± 0.6 nb, where the first uncertainty is statistical and the second is systematic. 
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The measurements of differential cross sections for heavy-quark production in high-energy 
hadronic interactions are critical input for the underlying next-to-leading order (NLO) Quan- 
tum Chromodynamics (QCD) calculations [1]. While progress has been achieved in the under- 
standing of heavy-quark production at Tevatron energies Il2l- [T0| , large theoretical uncertainties 
remain due to the dependence on the renormalization and factorization scales. Measurements 
of b-hadron production at the higher energies provided by the LHC represent an important 
new test of theoretical approaches that aim to reduce the scale dependence of NLO QCD calcu- 
lations [TTTT fT2| - The Compact Muon Solenoid (CMS) experiment, that covers a rapidity range 
complementary to the specialised b-physics experiment LHCb [13], recently measured the cross 
sections for production of B + [14] and B° [15] in pp collisions at y^s = 7 TeV. This paper presents 
the first measurement of the production of Bg, with Bg decaying into ]/ip (p, and adds informa- 
tion to improve the understanding of b-quark production at this energy. Data and theoretical 
predictions are compared at the level of NLO heavy-quark production. 

The decay channel Bg — > ]/ip cp is of wide interest for the indirect search of physics beyond the 
standard model at the LHC. This decay proceeds via the b — > ccs transition that probes the CP- 
violating phase related to Bg-Bg mixing. The standard model predicts this phase to be close to 
zero [16] while new phenomena may alter the observed phase l[T7 [ . The differential production 
cross section analysis prototypes the signal extraction towards time-dependent measurements 
in this final state. 

A sample of exclusive Bg — > ]/ip cp decays, with J/tp — > and <p —¥ K + K , is reconstructed 

in the data collected in 2010 by the CMS experiment, corresponding to an integrated lumi- 
nosity of 39.6 ± 1.6 pb -1 . The differential production cross sections, da/ dp\ and da/ dy B , are 
determined as functions of the transverse momentum p B and rapidity |y B | of the reconstructed 
Bg candidate. Here, the rapidity y B is defined as \ In §3^, where E is the particle's energy 
and Pl is the particle's momentum along the counterclockwise beam direction. The differential 
cross sections are calculated from the measured signal yields (n S i g ), corrected for the overall 
efficiency (e), bin size (Ax, with x = p\, |y B |), and integrated luminosity (L), 



where B is the product of the branching fractions for the decays of the J/tp and (p mesons. In 
each bin the signal yield is extracted with an unbinned maximum likelihood fit to the ]/\p <* 



invariant mass and proper decay length ct of the Bg candidates. The factor of 2 in Eq. ([lj 
accounts for the choice of quoting the cross section for Bg production, while n S i g includes both 

Bg and B^ states. The size of the bins is chosen such that the statistical uncertainty on n S i g is 
comparable in each of them. 

A detailed description of the CMS detector can be found elsewhere |[T8"l The primary compo- 
nents used in this analysis are the silicon tracker and the muon systems. The tracker operates in 
a 3.8 T axial magnetic field generated by a superconducting solenoid having an internal diam- 
eter of 6 m. The tracker consists of three cylindrical layers of pixel detectors complemented by 
two disks in the forward and backward directions. The radial region between 20 and 116 cm is 
occupied by several layers of silicon strip detectors in barrel and disk configurations, ensuring 
at least nine hits in the pseudorapidity range \rj\ < 2.4, where rj = — In [tan (9/2)] and 9 is the 
polar angle of the track relative to the counterclockwise beam direction. An impact parameter 
resolution around 15 pm and a px resolution around 1.5% are achieved for charged particles 
with transverse momenta up to 100 GeV/c. Muons are identified in the range \r/\ < 2.4, with 
detection planes made of drift tubes, cathode strip chambers, and resistive plate chambers, 



da(pp -> B° -> ]/tf> <p) 



(1) 



dx 



2 ■ e ■ B ■ L ■ Ax 




2 



embedded in the steel return yoke. 

The first level of the CMS trigger system uses information from the crystal electromagnetic cal- 
orimeter, the brass/ scintillator hadron calorimeter, and the muon detectors to select the most 
interesting events in less than 1 fis. The high level trigger (HLT) employs software algorithms 
and a farm of commercial processors to further decrease the event rate using information from 
all detector subsystems. The events used in the measurement reported in this paper were col- 
lected with a trigger requiring the presence of two muons at the HLT, with no explicit momen- 
tum threshold. 

Reconstruction of Bg — > J/ip (p candidates begins by identifying ]/xp — > decays. The muon 

candidates must have one or more reconstructed segments in the muon system that match 
the extrapolated position of a track reconstructed in the tracker. Furthermore, the muons are 
required to lie within a kinematic acceptance region defined as: pj > 3.3 GeV/ c for |//'' | < 1.3; 
total momentum p?' > 2.9 GeV/c for 1.3 < |^| < 2.2; and p\ > 0.8 GeV/c for 2.2 < \tjV\ < 
2.4. Two oppositely charged muon candidates are paired and are required to originate from 
a common vertex using a Kalman vertex fit. The muon pair is required to have a transverse 
momentum pj > 0.5 GeV/ c and an invariant mass within 150 MeV/ c 2 of the world average ]/ip 
mass value ||T9| , which corresponds to more than three times the measured dimuon invariant 
mass resolution 1120 1 . 

Candidate <p mesons are reconstructed from pairs of oppositely charged tracks with pj > 
0.7 GeV/ c that are selected from a sample with the muon candidate tracks removed. The tracks 
are required to have at least five hits in the silicon tracker detectors, and a track x 2 per degree 
of freedom less than five. Each track is assumed to be a kaon and the invariant mass of a track 
pair has to be within 10 MeV/ c 2 of the world average <£-meson mass |T9fl. 

The candidates are formed by combining a ]/xp with a (p candidate. The two muons and 
the two kaons are subjected to a combined vertex and kinematic fit |2T1| , where in addition the 
dimuon invariant mass is constrained to the nominal ]/\p mass. The selected candidates must 
have a resulting x 2 vertex probability greater than 2%, an invariant mass between 5.20 and 
5.65 GeV/c 2 , and be in the kinematic range 8 < pj < 50 GeV/c and |y B | < 2.4. For events 
with more than one candidate, the one with the highest vertex-fit probability is selected, which 
results in the correct choice 97% of the time, as determined from simulated signal events. 

The proper decay length of each selected Bg candidate is calculated using the formula ct = 
c(Mb/ Pj)L X y, where the transverse decay length L xy is the length of the vector s pointing from 
the primary vertex H22| to the secondary vertex projected onto the Bg transverse momentum: 
L X y = (s • p®)/pj, with Mb being the reconstructed mass of the Bg candidate. Candidate Bg 
mesons are accepted within the range —0.05 < ct < 0.35 cm. 

A total of 6200 events pass all the selection criteria. The efficiency of the Bg reconstruction is 
computed with a combination of techniques using the data and large samples of simulated 
signal events generated using pythia 6.422 l23l . The decays of unstable particles are described 
by the evtgen |l2"4"l simulation. Long-lived particles are then propagated through a detailed 
description of the CMS detector based on the GEANT4 [25J package. The trigger and muon- 
reconstruction efficiencies are obtained from a large sample of inclusive ]/xp — > decays in 

data using a (tag-and-probe) technique similar to that described in Ref. [20}. where one muon 
(the tag) is identified with stringent quality requirements, and the second muon (the probe) 
is identified using information either exclusively from the tracker (to measure the trigger and 
muon-identification efficiencies), or from the muon system (to measure the silicon tracking 
efficiency). The dimuon efficiencies are calculated as the product of the single-muon efficiencies 
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obtained with this method. Corrections to account for correlations between the two muons (1- 
3%) are obtained from simulation studies. The correction factors are determined in bins of 
single muon pj and t]V, and are applied independently to each muon from a Bg — > }/ip <p decay 
in the simulation to determine the total corrected efficiency. The probabilities for the muons to 
lie within the kinematic acceptance region and for the (p and Bg candidates to pass the selection 
requirements are determined from the simulated events. The efficiencies for hadronic track 
reconstruction 112611 and the vertex-quality requirement are found to be consistent between real 
data and simulated events within their uncertainties (up to 5%). The total efficiency of this 
selection, defined as the fraction of Bg — > ]/ip cp decays produced with 8 < pj < 50 GeV/ c and 
\y B | < 2.4 that pass all criteria, ranges from 1.3% for p\ « 8 GeV/c to 19.6% for p\ > 23 GeV/c. 

Backgrounds are dominated by prompt and non-prompt ]/xp production. The latter back- 
ground is dominated by B + and B° mesons that decay to a }/ip and a higher-mass K-meson state 
(such as the K + i). Such events tend to contribute to the low-mass side of the Mb mass distri- 
bution. Inspection of a large variety of potential background channels confirms that there is no 
single dominant component and that the channel B° -> J/xp K*(892) (with K*(892)° -> K+7T~), 
which a priori is kinematically similar to the signal decay and more abundantly produced, is 
strongly suppressed by the restriction on the K + K invariant mass. A study of the sidebands of 
the dimuon invariant mass distribution confirms that the contamination from events without a 
J/t/> decay to two muons is negligible after all selection criteria have been applied. 

The signal yields in each pj and |y B | bin, given in Table [lj are obtained using an unbinned 
extended maximum-likelihood fit to Mb and ct . The likelihood for event j is obtained by sum- 
ming the product of the yield n, and the probability density functions (PDF) Vi and Qi for 
each of the signal and background hypotheses i. Three individual components are considered: 
signal, non-prompt b — > ]/ip X, and prompt J/i/>. The extended likelihood function is then the 
product of likelihoods for each event j: 




The PDFs Vi and Q ; are parameterized separately for each fit component with shape param- 
eters a, for Mb and for ct. The yields n,- are then determined by minimizing the quantity 
— In C with respect to the signal yields and a subset of the PDF parameters |2"7|. Possible cor- 
relations between Mb and ct are found to be less than 2%. Therefore, they are assumed to have 
a negligible impact on the fit, and potential biases arising from this assumption are accounted 
for in the systematic uncertainty on the fitted signal yield as described below. 

The PDFs are constructed from basic analytical functions that satisfactorily describe the vari- 
able distributions from simulated events. Shape parameters are obtained from data when pos- 
sible. The Mb PDF is the sum of two Gaussian functions for the signal, a second-order polyno- 
mial for the non-prompt ]/\p that allows for possible curvature in the shape, and a first-order 
polynomial for prompt ]/ip. The resolution on Mb is approximately 20 MeV/c 2 near the Bg 
mass. 

For the signal, the ct PDF is a single exponential parameterized in terms of a proper decay 
length ct. It is convolved with a resolution function that is a combination of two Gaussian func- 
tions to account for a dominant core and small outlier distribution; the core fraction is varied in 
the fit and found to be consistently larger than 95%. The ct distribution for the non-prompt J/i/? 
background is described by a sum of two exponentials, with effective lifetimes that are allowed 
to be different. The "long-lifetime exponential" corresponds to decays of b-hadrons to a ]/xp 
plus some charged particles that survive the (p selection, while the "short-lifetime exponential" 
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Table 1: Signal yield n S i g/ efficiency e (%), and measured differential cross sections da/dp® and 
da/dy B , compared to the MC@NLO and pythia predictions, in different p B and \y B \ intervals. 
The uncertainties on the measured cross sections are statistical and systematic, respectively, 
excluding the common luminosity of 4% and the 1.4% from the ]/rp and (p branching fractions. 



p? (GeV/c) 


^sig 


e (%) 


da/dp B (nb/ GeV/c) 










Data 


MC@NLO 


PYTHIA 


8-12 


138 ± 16 


1.28 ± 0.05 


1.172 ±0.136 ±0.113 


0.719 


1.513 


12-16 


176 ± 17 


5.26 ± 0.23 


0.364 ± 0.035 ± 0.034 


0.240 


0.515 


16-23 


162 ± 16 


11.9 ± 0.6 


0.085 ± 0.008 ± 0.008 


0.074 


0.144 


23-50 


86 ± 11 


19.6 ± 1.1 


0.007 ±0.001 ±0.001 


0.008 


0.010 


|y B l 


n sig 


£(%) 


da/dy B 


(nb) 










Data 


MC@NLO 


PYTHIA 


0.00-0.80 


151 ± 15 


2.75 ± 0.09 


1.484 ±0.147 ±0.148 


1.040 


2.281 


0.80-1.40 


144 ± 15 


4.65 ± 0.18 


1.123 ±0.117 ±0.102 


1.023 


2.051 


1.40-1.70 


129 ± 15 


5.68 ± 0.31 


1.634 ±0.190 ±0.160 


0.929 


1.833 


1.70-2.40 


139 ± 17 


3.26 ± 0.20 


1.316 ±0.161 ±0.139 


0.801 


1.559 



accounts for events where the muons from the J/tp decay are wrongly combined with hadron 
tracks originating from the pp collision point. The exponential functions are convolved with a 
resolution function with the same parameters as the signal. For the prompt ]/ip component the 
pure resolution function is used. The core resolution in ct is measured in data to be 45 fim. 

All background shapes are obtained directly from data, while the signal shape in Mb is taken 
from a fit to reconstructed signal events from the simulation. The effective lifetime and resolu- 
tion function parameters for prompt and non-prompt backgrounds are extracted, using the full 
data sample irrespectively of p B and \y B \, from regions in Mb that are separated by more than 
four times the width of the observed B° signal from the mean Bg peak position (Mb sidebands): 
5.20 < M B < 5.29 GeV/c 2 and 5.45 < M B < 5.65 GeV/c 2 . A comparison of the PDF shapes for 
the different sideband regions in simulated events confirms that their average over the signal- 
free regions is a good representation of the background in the signal region. With the lifetimes 
for signal and non-prompt background fixed from this first step, the resolution function pa- 
rameters are then determined separately in each p B and \y B \ bin, from the Mb sidebands. The 
signal and background yields in each p B and \y B \ bin are determined in a final iteration, using 
the full Mb range, with all parameters floating except the background lifetimes and the lifetime 
resolution functions, which are fixed to the results of the fit to the Mb sidebands. It has been 
verified that leaving all parameters floating changes the signal yield by an amount smaller than 
the systematic uncertainty assigned to the fit procedure. 

Many detailed studies have been conducted to validate the accuracy and robustness of the fit 
procedure. A large number of pseudo-experiments were performed, each corresponding to the 
yields observed in each p B and \y B \ bin for a data sample corresponding to an integrated lu- 
minosity of 40 pb -1 , where signal and background events were generated randomly from the 
PDFs in each bin. The fit yields were found to be unbiased and their uncertainties estimated 
properly. The effects of residual correlations between Mb and ct were studied by mixing fully 
simulated signal and background events to produce pseudo-experiments. The observed devi- 
ations between the fitted and generated yields (1-2%) are taken as the systematic uncertainty 
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due to potential biases in the fit method. 

Figure [l] shows the fit projections for Mb and ct from the inclusive sample with 8 < p\ < 
50 GeV/c and |y B | < 2.4. When plotting Mb, the selection ct > 0.01 cm is applied for better 
visibility of the individual contributions. The number of signal events in the entire data sample 
is 549 ± 32, where the uncertainty is statistical only. The obtained proper decay length of the 
signal, ct = 478 ± 26 p.m, is within 1.4 standard deviations of the world average value IH9| , 
even though this analysis was not optimized for lifetime measurements. 

Table [l] summarizes the fitted signal yield in each bin of p B and | y B | . The differential cross 
section is calculated according to Eq. ([I]), using the product of the branching fractions B(]/ip — > 
p+}i-) = (5.93 ± 0.06) x 10~ 2 and B{(p -> K+K ) = (48.9 ± 0.5) x 10~ 2 HU. AU efficiencies 
are calculated separately in each bin, and account for bin-to-bin migrations (less than 1%) due 
to the finite resolution of the measured momentum and rapidity. 

The cross section measurement is affected by several sources of systematic uncertainty arising 
from uncertainties on the fit, efficiencies, branching fractions, and integrated luminosity. In 
every bin the total uncertainty is about 11%. Uncertainties on the muon efficiencies from the 
trigger, identification, and tracking are determined directly from data (3-5%). The uncertainty 
of the method employed to measure the efficiency in the data has been estimated from a large 
sample of full-detector simulated events (1-3%). The tracking efficiency for the charged kaons 
has been shown to be consistent with simulation. A conservative uncertainty of at most 9% 
in each bin has been assigned for the hadronic track reconstruction (adding linearly the un- 
certainties on the two kaon tracks II26II ), which includes the uncertainty due to misalignment 
of the silicon detectors. The uncertainty on the fit procedure arising from potential biases and 
imperfect knowledge of the PDF parameters is estimated by varying the parameters by one 
standard deviation (2-4%). The contribution related to the Bg momentum spectrum (1-3%) 
is evaluated by reweighting the shape of the p B distribution generated with pythia to match 
the spectrum predicted by MC@NLO [28J. An uncertainty of 1% is assigned to the variation 
of the selection criteria applied to the vertex fit probability, the transverse momentum of the 
kaons, the Bg transverse momentum, and the K + K invariant mass window. An uncertainty 
is added to account for the limited number of simulated events (at most 3% in the highest p B 
bin). The total uncorrelated systematic uncertainty on the cross-section measurement is com- 
puted in each bin as the sum in quadrature of the individual uncertainties, and is summarized 
in Table [l] In addition, there are common uncertainties of 4% from the integrated luminosity 
measurement [29 j and 1.4% from the ]/xp and (p branching fractions. As the reported result 
is a measurement of the Bg cross section times the Bg — > ]/ip (p branching fraction, the 30% 
uncertainty on the Bg — > ]/tp cp branching fraction [19J is not included in the result. 

The differential cross sections times branching fraction as functions of p B and |i/ B | are listed in 
Table [ljand plotted in Fig. [2} together with predictions from MC@NLO and PYTHIA. The predic- 
tions of MC@NLO use the renormalization and factorization scales ft = y m^c 4 + p T c 2 , where 

pi is the transverse momentum of the b quark, a b-quark mass of = 4.75 GeV/ c 2 , and 
the CTEQ6M parton distribution functions II30I . The uncertainty on the MC@NLO cross section 
is obtained simultaneously varying the renormalization and factorization scales by factors of 
two, varying by ±0.25 GeV/ c 2 , and using the CTEQ6.6 parton distribution function set. The 
prediction of pythia uses the CTEQ6L1 parton distribution functions 1130 1 , a b-quark mass of 
4.8 GeV/ c 2 , and the 72 tune ||3T| to simulate the underlying event. The total integrated Bg cross 
section times Bg — > ]/tp (p branching fraction for the range 8 < p T < 50 GeV/c and |y B | < 2.4 
is measured to be 6.9 ± 0.6 ± 0.6 nb, where the first uncertainty is statistical and the second is 
systematic. The statistical and systematic uncertainties are derived from the bin-by-bin uncer- 
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Figure 1: Projections of the fit results in M B (a) and ct (b) for 8 < pj < 50 GeV/c and \y B \ < 2.4. 
The curves in each plot are: the sum of all contributions (solid line); signal (dashed); prompt ]/xp 
(dotted); and non-prompt ]/xp (dot-dashed). For better visibility of the individual contributions, 
plot (a) includes the requirement ct > 0.01 cm. 
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Figure 2: Measured differential cross sections da I dpj (a) and da/ dy B (b) compared with the- 
oretical predictions. The (yellow) band represents the sum in quadrature of statistical and 
systematic uncertainties. The dotted (red) line is the PYTHIA prediction; the solid and dashed 
(blue) lines are the MC@NLO prediction and its uncertainty respectively. The common uncer- 
tainties of 4% on the data points, due to the integrated luminosity and of 30% on the theory 
curves, due to the Bg — > ]/xp (p branching fraction, are not shown. 
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tainties and propagated through the sum. The result lies between the theoretical predictions 
of MC@NLO (4.6+^7 ± 1.4 nb) and PYTHIA (9.4 ± 2.8 nb), where the last uncertainty is from 
the Bg — > ]/tp (p branching fraction fl9\. Also the previous CMS cross-section measurements 
of B + [14] and B° [15J production in pp collisions at \fs = 7 TeV, gave values between the two 
theory predictions, indicating internal consistency amongst the three different B-meson results. 

In summary the first measurements of the Bg differential cross sections da/ dpj and da/ dy B , in 
the decay channel B° — > ]/ip cp and in pp collisions at y/s = 7 TeV, have been presented. The 
results cover the kinematical window |y B | < 2.4 and 8 < pj < 50 GeV/c. They add com- 
plementary information to previous results in moving towards a comprehensive description of 
b-hadron production at y/s = 7 TeV. 
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